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ABSTRACT 

Using a sample o f 425 nearby Brightest Cluster Galaxies (BCGs) from 
Ivon der Linden et al.l (|2007ll , we study the relationship between their internal proper¬ 
ties (stellar masses, structural parameters and morphologies) and their environment. 
More massive BCGs tend to inhabit denser regions and more massive clusters than 
lower mass BCGs. Furthermore, cDs, which are BGGs with particularly extended 
envelopes, seem to prefer marginally denser regions and tend to be hosted by more 
massive halos than elliptical BCGs. cD and elliptical BCGs show parallel positive 
correlations between their stellar masses and environmental densities. However, at a 
fixed environmental density, cDs are, on average, ~ 40% more massive. Our results, 
together with the findings of previous studies, suggest an evolutionary link between 
elliptical and cD BCGs. We suggest that most present-day cDs started their life as 
ellipticals, which subsequently grew in stellar mass and size due to mergers. In this 
process, the cD envelope developed. The large scatter in the stellar masses and sizes 
of the cDs reflects their different merger histories. The growth of the BGGs in mass 
and size seems to be linked to the hierarchical growth of the structures they inhabit: 
as the groups and clusters became denser and more massive, the BCGs at their cen¬ 
tres also grew. This process is nearing completion since the majority (~ 60%) of the 
BCGs in the local Universe have cD morphology. However, the presence of galaxies 
with intermediate morphological classes (between ellipticals and cDs) suggests that 
the growth and morphological transformation of some BGGs is still ongoing. 

Key -words: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — 
galaxies: evolution — galaxies: formation 


1 INTRODUCTION 


The brightest cluster galaxies (BCGs) are the most lu¬ 
minous and massive galaxies in the universe. They are 
found at the centres of galaxy clusters and groups, 
and e xhibit many unique properti e s (see, e.g 


19871: Kormendv fc Diorrovskil Il989l : I Jordan et ahl l20d 


Ton^ 
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von der Linden et al. 1120071). Their origin and evolution is 
intimately linked with the evolution of their host clusters, 
and therefore can provide direct information on the for- 
mation and history of large-scale structures in Universe 
llConrov et ahllioo? !. 

Many scenarios have been proposed to explain the 
formation and e volution of BCGs. One of them is galac- 
tic cannibalism ilWhitel 1 19761 : lOstriker fc Hausrn^ ]l977l : 


* E-mail: : ppxdzl@nottinghain.ac.uk 
f E-mail: : alfonso.aragon@nottinghain.ac.uk 
X E-mail: : christopher.conselice@nottinghain.ac.uk 


iGariio et al.lll99^ . where BCGs were formed as a result of 
hierarchical mergers of smaller galaxies. O ther hypotheses 
include tidal stripping from cluster galaxies (jRichst ondl 19 76l : 
Im errittllig^ l. and star formation in the cluster core, w here 
BCGs are formed through cooling flows (lFabiarJll994h . Re¬ 
cently, numerical simulations and semi-analytic models sug¬ 
gest a two-phase process for BCGs formation. In these mod¬ 
els, the stellar component of BCGs was initially formed 
through the collapse of cooling gas or gas-rich mergers at 
high redshifts; subsequently, BCGs continued to grow sub- 
stantially by dissipationle s s processes such as dry mergers 
llDe Lucia fc Blaizotl l2007l : iNaab et"^ l2009l : iLaporte et al.l 
l2012h . This inside-out formation scenario is broadly consis¬ 
tent with observations, avoiding the need for cooling flows 
to provide the cold gas that would be necessary if BCGs had 
formed at later times. It also overcomes the problem caused 
by the merger rate in clusters being too low due to the high 
velocity dispersion in dynamically r elaxed clusters. However, 
some studies such as lAscaso et al.l 1 I 2 OIIII claimed that feed- 
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back rather than merging processes are the main mechanism 
affecting the evolution of the BCGs to the present epoch, 
ending the star formation within these systems. Therefore, 
many important details in the processes governing BCG for¬ 
mation and evolution are still unclear and deserve further 
investigation. 

Since BCGs posses singular properties (e.g., distinct 
structures and morphologies, and very high stellar masses) 
and reside in special environments (the core of groups and 
clusters), studying the relationship between their properties 
and their environments will help to constrain the theories of 
BCG formation and evolution and tell us whether the intrin¬ 
sic properties of BCGs or the environment play a dominant 
role in their history. In this context, it is important to bear 
in mind that, while both the location of BCGs at the bottom 
of the potential wells of clusters and their dominance at the 
massive end of the galaxy luminosity function may influence 
their properties, it is nonetheless very difficult to disentan¬ 
gle these two influences since it is hard to find equally mas¬ 
sive non-BCGs for comparison. Therefore, when comparing 
BCGs and non-BCGs, differences in the mass range spanned 
by the samples may bias the results. 

One key observational property of BCGs is that many 
of them sh ow unique morph ologies. The vast majority (but 
not all, see lZhao et al.ll201!Th BCGs are early-type galaxies. 
Mo st BCGs are classified as either elliptical or cD galax¬ 
ies dhauer fc PostmarJllQ^ : iFasano et al.ll201^ : IZhao et ^ 
120151 1. The defining characteristic separating these two mor¬ 
phological types is the presence of an extended, low-surface- 
brig htness stellar envelope in cDs that is absent in ellipticals 
fe.g. IPressleilllOsd : lOegerle fc HilllfiooH l. Since cDs are not 
found outside the BCG galaxy population, it is very im¬ 
portant to consider this unique galaxy class when studying 
BCGs. We will therefore use morphology as one of the main 
observables in this paper, focusing on the different proper¬ 
ties of elliptical and cD BCGs. 

Many previous observational works usually study the 
BCG population as a whole, and compare it with the popula¬ 
tion o f elliptical gala x ies th at are not BCGs ([Bernardi ^ aP 


20071: lLaueretal.1 l2007l: | von der Linden et al.l l2007l : 


Liu et all 2008tl . However, there has been some recent work 


exploring the structural differences between cluster ellipti¬ 
cals a nd BCGs with different morphologies. iFasano et al.l 
ll2010ll found that, while non-BCG cluster ellipticals 
generally have triaxial shape with a weak preference for 
prolateness, BCGs are also triaxial but with a much higher 
tendency towards prolateness. Such a strong prolateness 
appears entirely due to the fact that cDs dominate the 
BCG population. In fact, while the shape of elliptical BCGs 
does not differ from other cluster ellipticals, cDs tend 
to have prolate shapes. Furthermore, they suggest that 
the prolateness of the cDs could reflect the s hape of the 
associ ated dark matter haloes. More recently, IZhao et al.l 
(|2niRh have studied in detail the morphology and structure 
of BCGs, demonstrating that the morphological distinction 
between ellipticals and cDs is accompanied by quantitative 
structural differences. cD BCGs generally have much larger 
sizes and their light profiles cannot be modelled accurately 
using single Sersic functions. Conversely, elliptical BCGs 
are smaller and single Sersic profiles provide better fits to 
their surface brightness distributions. These differences in 
morphology and structure suggest that cD and elliptical 


BCGs have followed different evolutionary paths. We 
investigate these possible scenarios in this paper. 

There has been a significant amount of work address- 
ing the fo r matio n and evolution of BCGs. For example, 
I Guo et ^ (120091 ') studied how the structural parameters of 
central cluster galaxies correlate with their stellar masses 
and their host dark matter (DM) halo mass. They found that 
stellar mass is the dominant property dictating the shape 
and size of these galaxies, and suggest that the DM hal o 
mass does not pla y a ver y sign ificant role. HoggetaL JgOOdll , 
iKauffmann et ahl ll2004) and Ivan der Wei et al.l 1 200^ also 
reached similar co nclusions. In contrast, other studies (e.g., 
lAscaso et al. liml) claimed that there is a significant corre¬ 
lation betwee n the cluster mass and the p r opert ies of BCGs. 
Furthermore, pToymassian fc Andernachl (l2012l 'l added the 
cluster richness to the halo/cluster mass as another envi¬ 
ronmental indicator. They found that the absolute A-band 
luminosity of cD galaxies (a good proxy for stellar mass) 
strongly depends on the cluster richness, but less strongly on 
the cluster velocity dispersion (a proxy for DM halo mass). 
Therefore, since the effects of the halo mass and the cluster 
richness could be different, it is necessary to take them into 
account as separate environmental parameters when study¬ 
ing BCG evolution. 

Many other recent papers have studied the prop¬ 
erties of BCGs in relation to other early-type galaxies, 
providing important c lues to how they form an d evolve . 
Some examples include |Shanka r et al.l (I2 013l.l2014af9. [2015l f : 
iHuertas-Companv et al.! ^ 2013a Bemardif i 2009 ). For the 
sake of brevity, we will not describe their findings here but 
we will mention them in the following discussion when rele¬ 
vant. 

In this paper we use a well-defined local sample of 625 
BCGs from Ivon der Linden et ^ (l2007l . hereafter L07) and 
carry out a comprehensive and systematic statistical study 
on the correlation between BCGs intrinsic properties (struc¬ 
ture, morphology and stellar mass) and their environment. 
We consider two environmental measures, a global one (the 
DM cluster halo mass, characterised by its velocity disper¬ 
sion) and a local one (the galaxy density). In doing so we will 
obtain very valuable additional information on how BCGs 
form and evolve. 

The galaxy groups and clusters these BCGs inhabit 
span a very broad range of total masses, from ~ Mq 
to ~ Since there is no clear boundary separat¬ 

ing “clusters” from “groups” (although 10^^ Mq could be 
taken as the transition mass), we will study group and clus¬ 
ter BCGs together. We will explore how the masses of the 
parent groups/clusters affect the properties and evolution of 
the BCGs. 

The paper is organized as follows. In (J21 we introduce 
the BCG sample, and describe the observables we will use 
(morphologies, structural parameters, stellar masses, envi¬ 
ronmental densities, and DM halo virial masses). In 321 we 
show how the structural parameters of the BCGs relate to 
their stellar masses, and their global and local environment, 
and discuss the implications of the correlations we find on 
the formation of the BCG population. In 31] we go one step 
further and bring the galaxy morphologies into the general 
picture to learn about the distinct evolutionary history of cD 
and elliptical BCGs. We summarise our main conclusions in 
31 Throughout this paper we have adopted the ACDM cos- 
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mology with = 0.3, JIa = 0.7, and Hq — 70 km s ^ 
Mpc"^ 


2 BCG SAMPLE AND PROPERTIES 

The parent BCG sample we use in this paper comes from the 
catalogue published by L07. The groups and clusters that 
host these BC Gs are contained in the SDSS-based C4 clus¬ 
ter catalogue llMiller et al.ll2005l l. a widely-used and well- 
defined sample whose reliability has been thoroughly tested 
by simulations. Based on the C4 sample, L07 developed 
an improved algorithm to identify the BGG in each clus¬ 
ter and published a catalogue containing 625 BCGs resid- 
ing in galajcy groups and c lusters at 0.02 ^ a ^ 0.10. See 
Ivon der Linden et al. Idio^ for a detailed discussion on the 
BCG identification method. _ 

In our previous paper llZhao et al.ll201^ we published 
visual morphologies for these 625 BCGs. The BCGs were 
classified into three main types: 414 cDs, including pure cD 
(356), cD/E (53) and cD/SO (5); 155 ellipticals, including 
pure E (80), E/cD (72), and E/SO (3); 46 disk galaxies, con¬ 
taining spirals (24) and SOs (22). There are also 10 BGGs un¬ 
dergoing major mergers. We used intermediate classes such 
as cD/E (probably a cD, but could be E) and E/cD (prob¬ 
ably E, but could be cD) to account for the uncertainty in¬ 
herent in the visual classification. Separating cD BCGs and 
non-cD elliptical BCGs is a very hard problem since there is 
no sha rp morphological distinction betwe en these two classes 
(e.g., IPatel etld] 20061 : llTu et akllioosh . Detecting the ex¬ 
tended stellar envelope that characterises cD galaxies de¬ 
pends not only on its dominance, but also on the quality and 
depth of the images, and on the details of the classification 
method(s) employed. Since the SDSS imaging data that we 
use is of uniform quality, the c lassifi cation is internally con¬ 
sistent. Moreover. I Zhao et al. I diolsl) demonstrated that the 
morphological type of a galaxy is very tightly related to its 
structural parameters, and devised a quantitative method to 
separate cDs from elliptical BCGs in a robust manner that 
agrees very well with the visual classification. We are there¬ 
fore confident that the morphological information that we 
use here is reliable and self-consistent within our sample. Ob¬ 
viously, a degree of caution would be necessary when com¬ 
paring our morphologies with those of galaxies from other 
samples since the quality of the images and the classifica tion 
criteria may be different. However, in IZhao et al.l ll2015l l we 
carried out imaging simulations to show that the outer en¬ 
velopes of cD galaxies would have been detected if present, 
and thus that there is no bias in our morphological distinc¬ 
tion between E and cD. 

In this paper we w i ll call “cD BCGs” the 414 galaxies 
classified bv lzhao et al.l (l2015l i as cD, cD/E and cD/SO, and 
“elliptical BCGs” the 155 galaxies classified as E, E/cD and 
E/SO. We will also include in our study the 46 disk BCGs 
(spirals and SOs), but not the 10 major mergers. This sample 
therefore contains 615 BCGs. 

The structural properties (Sersic index n and effective 
radius i?e)EI that we use in this paper were also published 

^ Strictly speaking, Re is the effective semi-major axis of the 
single Sersic model fit 


bv izhao et al.l (l2015ll . These were derived from SDSS DR7 
r-band images using two-dimensional single ISersi 2 (Il963ll 
model fits to the galaxi es’ light profiles. The fits were car¬ 
ried out with GALFIT jPen g et al.ll20o3 l using the GALA¬ 
PAGOS llBarden et al.ll2012l i pindine. The method simul¬ 
taneously fits the target galaxy and its near neighbours, 
yielding more accurate fits and improved sky subtraction. 
The imaging data reaches a surface brightness limit of 27 
mag/arcsec^, and are therefore deep enough to study the 
faint extended envelopes present in cD BCGs. The detailed 
description of the fitting pro cedure and structu ral parame¬ 
ter estimation can be found in IZhao et al.l (l2015li . The values 
of Re and n that we obtained are broa dly compatible with 
the ones published bv iGuo et al.l (120091 '). However, there are 
some relatively minor systematic differences due to the im- 
pr ovements in the sk y subtraction procedure implemented 
by I Zhao et al.l (l2015li . A direct comparison is presented in 
Appendix \M 

The stellar masses we use come from “The MPA- 
JHU DR7 release of spectrum measurements” (see 
www.mpa-garching.mpg.de/SDSS/DR7/j3- Hereafter we call 
these “MPA-JHU masses”. These stellar masses are ob¬ 
tained via spectral energy distrib ution (SEP) fit s to the 
DR7 photometric data using a iKrou 3 (l200lf ) Initial 
Mass F unction. Although the m eth od is not iden t ical t o 
that of iKauffmann et al.l (l2003l l or iGallazzi et al.l ll2005h . 
who use spectroscopic information, the resulting masses 
agree very well with only a few minor offsets. A de¬ 
tailed discussion and comparison of the methods can be 
found in www.mpa-garching.mpg.de/SDSS/DR7/mass_comp. 
The number of BCGs in our sample which have MPA-JHU 
stellar mass information is 591, i.e., 96%. The very small mi¬ 
nority of galaxies without stellar masses include 20 galaxies 
for which no spectroscopic redshift is available (essential to 
determine accurate distances) and 4 for which the MPA- 
JHU catalogue fails to provide a value for the mass, pre¬ 
sumable because the SED fitting method does not yield a 
reliable solution. Since only 4% of the galaxies in the par¬ 
ent sample do not have stellar masses, we do not expect 
them to have any significant influence in our results. At this 
stage, and in order to ensure we have a stellar-mass-selected 
sample, we impose a minimum mass of 3 x Mq, which 
reduces the sample to 535 BCGs. This limit also eliminates 
a few galaxies whose stellar masses, structural parameters 
and morphologies have larger uncertainties due to their faint 
magnitudes. 

These MPA-JHU stellar masses are derived from Pet¬ 
rosian magnitudes and are therefore not dependent on the 
fitting parameters that we obtain. This is important since it 
allows us to look for independent correlatio ns between stellar 
mass and the fit parameters. Alternatively. f^uo et al.l (l2009li 
estimated stellar masses using photometric fluxes derived 
from their light profile model fits. Such a method results in 
model-dependent stellar masses, which may produce spuri¬ 
ous correlations between the masses and the model parame¬ 
ters. We will discuss this in more detail in section [S] and we 
will argue that for our study the MPA-JHU Petrosian-based 
stellar masses should be preferred. 


^ In this paper we use their updated stellar masses from 
http: //home. strw. leldenuniv. nl/~ j arle/SDSS/ 
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The final key ingredients in our study are quantitative 
measurements of the environments where the BCGs reside. 
We will use two distinct descriptions of the environment, 
global and local. The “global environment” is governed by 
the properties of the cluster/group that contains the BCG, 
and in particular its total mass (including the dark-matter 
halo). We use the velocity dispersion of the cluster (cr 2 oo) 
published by L07 t o estimate the halo virial mass M 200 using 
the Equation 10 of iFinn et al.l (l2005l l , which is 


A /200 = 


( loooZ,-. ' 


The group and cluster sample studied here covers a broad 
range of masses, from M 200 ~ 10^^ A/q to A /200 ~ 10^® A/©, 
peaking at A /200 ~ 10^"^ A/© (see Fig. [6|. 

To characterise the “local environment” we use the en¬ 
vironm ental luminosity density introduced bv lTempel et al.l 
(l2012t) . This is a good proxy for t he environmental stellar 
mass density, which, as argued bv IWolf et al.l ll2009tl . is a 
better and more robust measurement of the environment 
than galaxy number density. The main advantages of using 
stellar mass (or luminosity) density over galaxy number den¬ 
sity are twofold. First, the environmental luminosity/mass 
density does not depend strongly on the exact details of 
the galaxy sample used to define it, such as the magnitude 
limit, provided that it reaches significantly fainter than the 
“knee” of the luminosity function. And second, it represents 
better the strength of the interactions that a galaxy may 
experience from its neighbours: it is not the same to be sur¬ 
rounded by N faint low-mass galax ies than by N bright 
high-mass ones. iTempel et al.l (120121 ') determined these en¬ 
vironmental densities using SDSS r-band luminosities with 
a smoothing scale of lh“^Mpc. The total number of BCGs 
in our mass-limited sample for which we have both stellar 
masses and environmental densities is 425. The galaxies for 
which environmental densities are not available are outside 
the foo tprint of the contigu ous sky region covered by the 
work of iTempel et al.l (l2012l l. and therefore there is no rea¬ 
son to believe that their exclusion from our analysis will 
bias our conclusions. The BCG sample covers one order of 
magnitude in environmental density (see Fig. 1^. 

In what follows, we will consider the sample comprising 
the 425 A/* > 3 x M© BCGs with cD (275), elliptical 
(116), SO (15) and spiral (19) morphologies for which we 
have obtained stellar masses, cluster masses and environ¬ 
mental densities. 


3 CORRELATIONS BETWEEN BCG 
PROPERTIES 

In this section we analyse the correlations (or lack thereof) 
between the structural parameters, masses and environ¬ 
ments (global and local) of the BCG population as a whole 
and discuss their implications. In section |4] we will include 
morphology as an additional key property. 
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Figure 1. Comparison between the stellar masses and the struc¬ 
tural parameters of the BCGs in our sample. Upper panel: Sersic- 
index n vs. MPA-JHU stellar mass M*. Lower panel: effective 
radius Re vs. M*. Red plus signs, green crosses, magenta open 
squares and blue open triangles correspond to cD, elliptical, SO 
and spiral BCGs, respectively. Black solid squares with error bars 
in upper panel show the median and the 84 and 16 percentiles 
(~ Icr) of each parameter in 0.15 dex logM* bins for the com¬ 
bined cD and elliptical BCGs. Red dots and green diamonds with 
error bars in lower panel are for cD and elliptical BCGs, respec¬ 
tively. Bins with fewer than 20 galaxies are excluded due to their 
large statistical uncertainties. The black solid line in the lower 
panel corresponds to the best-fit relation for the normal (non- 
BCG ) early-type galax y population, defined to have n > 2.5, 
from IShen et al] J2003h . The dashed lines correspond to the Icr 
scatter in this relation. 


3.1 Stellar Masses and Structural Parameters 

First we explore the relation between the BCGs structural 
parameters (Sersic index n and effective radius Re) and their 
stellar mass M*. In the top panel of Fig. [J we investigate 
whether there is a statistical correlation between the galax¬ 
ies’ profile shape, characterised by n, and their stellar mass. 
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To guide the eye, we have binned the data in stellar mass 
bins 0.15 dex wide. The black squares with error bars show 
the median and the 84 and 16 percentiles Icr) of the n 
distributions for each mass bin, considering only the BCGs 
with cD and elliptical morphologies. In order to avoid large 
statistical uncertainties, we exclude bins with fewer than 20 
galaxies. 

We find no correlation between n and M* for these 
galaxies. The median n for the elliptical and cD BCGs is 
6.02, which indicates that, on average, these galaxies have 
both centrally-concentrated light profiles and extended en- 
velo pes, as expected for a population dominated by cDs 
(seelZhno_et_^ |20ld, and references therein). Interestingly, 
as IZhao et al.l ( 2015t ) pointed out, there is little separation 
between the n distributions of cD and elliptical galaxies. 
A Kolmogorov-Smirnov test indicates that the difference is 
only significant at the 2a level. The median Se rsic index n 
is 6.12jlj g 3 for cDs and 5.86lj'42 for ellipticals dZhao et al.l 
[20l5lf| . The slightly larger median n value of the cD galax¬ 
ies is driven by their extended envelope. As expected, disk 
BCGs (spirals and SOs) have significantly lower n values 
(2.91 and 3.88 respectively). 

The lack of correlation between n and for the 
BCGs in our sample contrasts with the findings of iGuo et al.l 
(l2009li . who claimed a clear positive correlation in the sense 
that more massive BCGs seem to have higher values of n. 
As we show in Appendix IA21 we believe this may be due 
to the fact that IGuo et al] ~( 2009l l estimated stellar masses 
from total luminosities derived from single Sersic model fits. 
These luminosities (and the derived stellar masses) depend 
on the value of n, and this dependency could drive an arti¬ 
ficial correlation. 

As an aside, we note that in the upper panel of Fig. [T] 
there is a small number of cD and elliptical BCGs whose n 
is quite large (n > 12). It is important to realise that for 
large n (n > 6 or so) very small changes in the light profile 
result in large changes in n, and thus all values of n above 
~ 6 correspond essentially to the same profile. Furthermore, 
a visual inspection of the fits and the residuals indicate that 
these large n objects are usually surrounded by multiple 
close bright companions (or, in a few cases, a bright nearby 
star). This makes the fits less reliable. Furthermore, some of 
these objects have double cores, and therefore a single Sersic 
profile is not a good model of their surface brightness distri¬ 
bution. In these cases, the derived model parameters should 
be taken with caution. Since the fraction of affected objects 
is quite small, they do not affect the statistical conclusions 
of this study. Removing them would have no significant sta¬ 
tistical effect, and they are therefore kept in our analysis 
for completeness. Another reason for this that the high n 
systems are distributed over all stellar masses, and not just 
found within the high or low stellar mass systems. 

We examine now the relationship between the effec¬ 
tive radius Re and the stellar mass of the BGGs shown 
in the lower panel of Fig. [T] For comparison, we show the 
relation found fo r normal non-BCG early-type galaxies by 
IShen et ^ (l2003l l selected from the SDSS survey as system 
with n > 2.5. The sizes and stellar masses published by 


® The errors quoted for median values correspond to the 84 and 
16 percentiles of the distributions (~ Icr scatter). 


IShen et aP ll2003l ') are directly comparable to the ones we 
use. Their effective radii are computed from single Sersic 
fits to SDSS images, like ours, an d their stellar ma s ses ar e 
also derived using the method of iKauffmann et al.l (l2003ll . 
Note that the lShen et al.l ll2003h sample is dominated by field 
galaxies, although we will see below that similar conclusions 
are obtained for cluster early-types. 

The effective radii of early-type BCGs is strongly cor¬ 
related with their stellar masses: on average, Re increases 
when M* increases, but the scatter is large (about ~ 0.3 dex, 
or a factor of 2 in Re at a given mass). In agreement 
with lSernar HI (l2009ll , we find that almost all the BCGs are 
above the average relation for non-BCG early types, and 
the slope is similar (within a large uncertainty). The scat¬ 
ter is also larger for the BCGs than for the other early-type 
galaxies. Notwithstanding this large scatter, the median ra¬ 
dius of BCGs is about twice as large as that of non-BCG 
early types of similar masses. This difference is largely due 
to the cD galaxies, which dominate the sample. As shown 
in the lower panel of Fig. [Jl when we analyse the proper¬ 
ties of BCGs separated by morphology, elliptical BCGs are, 
on average, significantly smaller than cDs. The minority of 
BCGs that have disk (spiral and SO) morphologies tend to 
populate the low end of the size distribution. 

Fig.[T] also shows that the BCGs in our sample span a 
very broad range of stellar masses (10^°'®-10^^ Mq)- This is 
mainly due to the fact that these BCGs are hosted by galaxy 
groups and clusters with very different masses (Fig. [6]), 
combined with the weak correlation between the galaxies’ 
stellar masse and M 200 (Fig. U). Nevertheless, it is clear 
that at all stellar masses BCGs have larger radii than non- 
BCG early-type galaxies . Thi s agrees with the fin dings of 
IValentinuzzi et ^ (l2010ll and IVulcani et akl (l2014ll for low- 
redshift BCG and non-BCG galaxies in the WINGS clusters 
(see their Fig. 11). Although a detailed quantitative com¬ 
parison is very difficult given the differences in methodology 
combined with the fact the the WINGS sample does not in¬ 
clude groups, it is reassuring to see that compatible results 
are obtained independently. Note also that the stellar masses 
of the WINGS BCGs are all in the range 
where most of our BCGs lie, but we also have BCGs with 
lower stellar masses since our sample includes both clusters 
and groups. 


3.2 Local environment: the effect of galaxy 
density 

We explore now the relationship between the local envi¬ 
ronment that BCGs inhabit and their intrinsic properties 
(structural parameters and stellar masses). As discussed in 
Section [H we use t he environmental luminosity density of 
iTempel et al.l (120121 ') to characterise the local environment. 
In the three panels of Fig. [2] we plot the Sersic index n, the 
effective radius Re, and the MPA-JHU stellar mass M, vs. 
this density. The left panel shows that there is no correlation 
between n and density (Pearson correlation coefficient 0.03). 
However, both Re and clearly correlate, on average, with 
density (correlation coefficients 0.32 and 0.49 respectively). 
Although there is significant scatter, larger and more mas¬ 
sive BGGs tend to inhabit in denser environments. 

It appears that local density correlates with both the 
size and the stellar mass of the early-type BCGs. However, 
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Figure 2. Relationship between environmental density and BCG properties. From left to right, these properties are the Sersic index n, 
the effective radius Re, and the stellar mass M«. Symbols as in Fig.[T] 




Figure 3. Left panel: Re vs. AL* for cD and elliptical BCGs binned by environmental density. Right panel: M* vs. environmental density 
binned by Re- The points correspond to the median for each bin and error bars indicate the 84 and 16 percentiles (~ Ict). Bins containing 
fewer than 5 galaxies have been excluded due to their large statistical uncertainties. The legend shows the different symbols corresponding 
to each bin. 


Fig. [T] shows that Re correlates with M*. It is therefore im¬ 
portant to ascertain which of these two parameters is the in¬ 
trinsic driver of the correlations with density. To do this, in 
the left panel of Fig. [3] we plot Re vs. M* binning the galax¬ 
ies by density. We only include cD and elliptical BCGs. For 
a given stellar mass, the median Re is the same for all den¬ 
sities. This suggests that density does not affect BCG size 
directly, but only through its dependence with stellar mass. 
In the right panel of this figure we show the M*-density 
relation again, but now binning the galaxies by radius. For 
galaxies of all sizes, there is a clear correlation between stel¬ 
lar mass and environment: more massive BCGs tend to in¬ 
habit denser regions, regardless of their radius. This implies 
that the stellar mass-density correlation is the more fun¬ 
damental one, and that the environment affects the BCG 
stellar mass more directly than their sizes. 

The fact that the mass-size relation for the general 
galaxy population does not depend signihcantly on en¬ 


vironment (at least at lo w redshift) has been found in 


several recent studies (e.g 

.. Shen et alJ 20031: Maltbv et alJ 

2 OIOI: Rettura et al. 2010l: 

1 1 

Huertas-Gomnanv et al.l 2013a bl: 

Poeeianti et al.||2013|). Our results reveal that this is also 


true for BCGs. 


3.3 Global environment: the effect of the cluster 
mass 

We now consider the effect of the global environment (char¬ 
acterised by the total mass of the host cluster M200', see 
Section [2J on the properties of the BCGs. Fig. [T] shows the 
relation of AL 200 with the Sersic index n, effective radius 
Re, stellar mass M* and environmental density (from left to 
right). 

The Sersic index does not show any dependence on 
the halo virial mass (Pearson correlation coefficient —0.04). 
Both effective radius and stellar mass show a small degree 
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Figure 4. Relationship between M 200 and other BCG properties. From left to right, these properties are the Sersic index n, the effective 
radius Re, the stellar mass M* and the environmental density. Symbols as in Fig. fTI 





log M ,/M g log g log density log density 

Figure 5. From left to right, the first panel shows Re vs. M« in M 200 bins; the second panel M* vs. M 200 in density bins; the third 
panel M* vs. density in M 200 bins; and the fourth panel M 200 vs. density in A4* bins. The points correspond to the median for each 
bin and error bars indicate the 84 and 16 percentiles (~ Itr). Bins containing fewer than 5 galaxies have been excluded due to their 
large statistical uncertainties. The legend shows the different symbols corresponding to each bin. Only cD and elliptical BCGs have been 
included. 


of correlation with M 200 , albeit with large scatter (correla¬ 
tion coefficients 0.26 and 0.17 respectively (Q. As before, we 
need to explore which of these two parameters is the driver 
of the observed correlations. The first panel of Fig. [S] shows 
that the stellar mass-size relation does not depend on the 
M 200 (global environme nt), in agreement with the findings 
of IShankar et al.l ll2014bl '). Since we also found in Section [3. 2 1 
that the size of BCGs is not directly affected by the local 
environment (or galaxy density) we conclude that any ap¬ 
parent environmental effect on Re is driven by the stellar 
mass-size relation combined with the environmental depen¬ 
dence (or dependencies) of stellar mass. 

We now consider the effect of environment on the BCGs’ 
stellar masses. Previous studies have found that the stellar 
masses of the BCGs correlate with the t otal mass (or veloc¬ 
ity dispersion) of the host cluster fe.g.. Iwhilev et ^l2008l : 
lAscaso et alj|201ll i. One complication that plagues all envi¬ 
ronmental studies is the fact that the two characterisations 
of the environment that we use (local and global) are, not 
surprisingly, correlated (see rightmost panel of Fig. HI), al- 

^ Note that Re and M« correlate more weakly with M 200 than 
with the environmental density (compare Figs. and |4)|. 


though not very tightly (correlation coefficient 0.33). How¬ 
ever, these two measures of environment are clearly not 
representing the same physical scales or the same range 
of ph ysical processes, and t heir evolution is largely decou¬ 
pled JPoggianti et al.l l2010l ~l . There is also clear evidence 
that local and global environment do not have the same ef - 
fect on galaxy evolution. For instance, IVulcani et al.l (l2012l l 
found that the local environment has a strong effect on the 
galaxies’ stellar mass function, while the same team showed 
that the global envi ronment has no (or much weaker) effect 
dVulcani et al.|[201^ 1. 

We find that the correlation between M* and environ¬ 
mental density (Fig. [2]right panel; Pearson correlation coef¬ 
ficient 0.49) is much stronger than the M,-M 2 oo one (Fig [4] 
third panel; correlation coefficient 0.17), suggesting that the 
main driver of these correlations is the local density. This is 
confirmed by Fig. [S] The second panel shows that at fixed 
density the correlation between and M 200 largely dis¬ 
appears, except, perhaps, for the two highest density bins, 
although the statistical uncertainties are large. However, the 
third panel indicates that at fixed M 200 the M,-density rela¬ 
tion is still present. The fourth panel shows that at fixed 
most of the M 2 oo-density correlation vanishes. We conclude 
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that the M,-environment correlations are really driven by 
the M*-density correlation, while the weaker M*-M 2 oo cor¬ 
relation is secondary, and it originates on the M 2 oo-density 
and M,-density correlations. 

It could be argued that the detected trend (more mas¬ 
sive BGGs live in denser, more massive halos) may be due, 
at least partially, to a pure statistical effect. If stellar masses 
are randomly drawn from the mass function of galaxies, mas¬ 
sive halos, which host a larger number of galaxies, have 


e.e:.. Tremaine & Richstonel 1977 

: Bhavsar & Barrow 19851: 

Lin et al. 2010l: iDobos & Gsabai 

I 2 OIII: Paraniaue & ShethI 

2 OI 2 I: lMorell2012|j. However, we 

argue that this statistical 


effect cannot be the main driver of the correlation we find. 
There is quite a lot of evidence indicating that the luminosity 
of cluster BGGs is inconsistent with just statistical sampling 
of the cluster galaxy luminosity function: BGGs are gener¬ 
ally too bright, and there is too large a gap betwee n the lu¬ 
minosity of the first and second brightest galaxies (ISandagj 
19761: Tremaine fc Richstond 19771: Bhavsar fc Barro^ 19851: 


Dobos fc Csabaill2011 ; Morel2012l : lllearin et al.l2013l . among 

others)0. 

If BGGs are not governed by the luminosity/mass func¬ 
tion of the rest of the cluster galaxies, the above statistical 
arguments do not apply. Things may be not so clear for the 
poorest groups, where the brightest galaxies seem to be com¬ 
patible with being statistically drawn from the bright end of 
the galaxy luminosity function, as argued by some of these 
authors. However, the correlation between BCG mass and 
environment appears stronger for more massive and denser 
clusters (see, e.g., rightmost panel of Fig[2]), where we ar¬ 
gue this statistical effect should not apply, and weaker for 
poorer groups, where the statistical bias should be strongest. 
If the main driver of the correlation were just the statistical 
sampling of the luminosity function, we would expect the 
correlation to be strongest where this effect is most impor¬ 
tant (low mass and less dense clusters and groups). Since the 
effect we find is strongest for high-mass and denser clusters, 
we conclude that the correlation cannot be primarily driven 
by sampling statistics. 


In summary, in this section we have found that BGGs fol¬ 
low a stellar mass-size relation that is independent of the 
environment, and that stellar mass is intrinsically correlated 
with the local environment (or environmental density). In 
Section |4] we will see how these correlations depend on the 
morphologies of the BGGs. 


4 EVOLUTIONARY HISTORY OF CD AND 
ELLIPTICAL BCGS 

In I Zhao et ^ (|2015h we found that the vast majority of 
BGGs (over 90%) have cD or elliptical morphologies, while 
only a small minority (-^ 7%) are disk galaxies (spirals and 


^ Note, how ever, that |PjnjMjjajje_fc_Sheth| ll2012l l disagree, but 
[m or i mH) and iHearin et~alTl[201.'ji hare argued against their 
results 


SOs), and the remaining few are major mergers. The mor¬ 
phology of these galaxies is clearly linked to their quanti¬ 
tative structural parameters. cDs are generally larger than 
ellipticals, and their light distributions deviate significantly 
more from Sersic profiles than those of ellipticals. With the 
additional information presented in this paper we will now 
explore how morphology and structure are linked to the stel¬ 
lar masses and environments of the BGGs. 

In Fig. [6] we present the distributions of the stel¬ 
lar masses, environmental densities and parent cluster to¬ 
tal masses (M 200 ) for cD and elliptical BGGs. The left 
panel clearly shows that cDs have, statistically, larger stellar 
masses than elliptical BGGs. The median stellar mass of the 
cDs is X IO^^Mq, ~ 50% larger than that of ellipti¬ 

cals (1.4lg g X IO^^Mq). a two-sample Kolmogorov-Smirnov 
test shows that this difference is significant at the 4.6 (t 
level. The disk galaxies (not shown in the figure for clarity) 
are even less massive: the median stellar mass for spirals and 
SOs is l.OtJ;^ X IO^Mq. 

With respect to environmental density (middle panel of 
Fig. [6l), cDs seem to prefer marginally denser regions (by 
~ 20% on average) than elliptical BGGs, although, statisti¬ 
cally, this difference is only significant at the ~ 2.4(j level. 
Disk galaxies tend to live in the regions with the small¬ 
est densities (a factor of ~ 2 smaller than cDs). Similarly 
(right panel of Fig. |6]), cDs appear to be hosted by more 
massive clusters/groups than ellipticals, but once again the 
difference (a factor of ~ 1.7 in median M 200 ) is only barely 
significant (~ 2a). 

These differences in the stellar masses and environments 
of BGGs with different morphologies suggest that their for¬ 
mation histories may be different. In Section[3]we found that 
there are intrinsic correlations between Rs and M*, and be¬ 
tween M* and the environmental density. By exploring the 
relationship between these properties and the galaxies’ mor¬ 
phologies we may be able to shed additional light on the issue 
of the formation and growth of BGGs. In Fig. [7] we show the 
Re-density relation (left panel) and the M*-density relation 
(right panel) for cD, elliptical, and disk BGGs. cD and el¬ 
liptical BGGs show parallel correlations, in the sense that 
larger and more massive galaxies tend to prefer denser en¬ 
vironments. However, at a fixed environmental density, cDs 
are, on average, a factor of ~ 2 larger and ~ 40% more mas¬ 
sive than elliptical BGGs. Disk galaxies tend to be smaller 
and less massive, but clear correlations are not seen, perhaps 
due to the small number statistics. This correlation is also 
seen when investigating the relation with the the total mass 
of the cluster. These correlations are futhermore certainly 
due to the fact that there is a different relation between the 
stellar mass and radius for ellipticals and cD. This effect is 
driven by the stellar mass being higher, which then increases 
the radius. 

Note that the observational results presented in this 
paper, including the differences found between cDs and el¬ 
liptical BGGs, do not depend on whether the morphological 
classihcation is done visually (as shown here) or automat- 
ically (based on structural parameter method described in 
IZhao et al. IlioTi). A parallel analysis using the automatic 
cD/elliptical classihcation yields entirely consistent results. 
We are therefore conhdent that our results are robust, and 
do not depend signihcantly on the details of the morpholog¬ 
ical classihcation. 
















































Evolution of the Brightest Cluster Galaxies 9 



Figure 6. Distribution of M*, environmental density and M 200 for the 275 cD (red solid) and 116 elliptical (green dashed) BCGs in our 
sample. The cr value in each panel indicates the significance (confidence level) of the observed differences between the cD and elliptical 
BCG parameter distributions. These are derived from two-sample Kolmogorov-Smirnov tests. Statistically, compared with elliptical 
BCGs, cD galaxies are more massive, tend reside in denser environments, and tend to be hosted by more massive dark matter halos. The 
median values of the different distributions are indicated by the vertical lines and adjacent numerical values. 



D8 ^1 ^4 ^7 

log density 


Figure 7. i?e-density and M*-density relations for BCGs with different morphologies. Red plus signs, green crosses, magenta open 
squares and blue open triangles correspond to cD, elliptical, SO and spiral BCGs, respectively (as in Fig. Red filled circles with error 
bars show the median and the 84 and 16 percentiles for cD galaxies. Green filled diamonds show the same properties for elliptical BCGs. 
The red and green lines show a linear fit for cD and ellipical BCGs Respectively. It is clear that at the same density, cD galaxies are 
statistically larger by factor of ~ 2 than elliptical BDGs. The stellar mass of cDs is larger by a factor of ~ 1.4 than that of ellipticals. 
Disk BCGs tend to be smaller and less massive. 



log density 


Our empirical results, together with the findings 
of previous works, suggest a possible sce nario linking 
the e vo lution of elliptic a l and cD BCGs. Whilev et alJ 
(|2008|), iBurke &; Collinsl (|2Q13l) . iBurke et al. ( 20151) and 
IZhang et al.l ( 20151) . among others, suggest that the stellar 
mass of BCGs has experienced some (but relatively mod¬ 
erate) growth in the last ~ 6-8 Gyrs. Although measuring 
BC G growth is notoriou sly difficult due to progenitor bias 
(see IShankar et alJ l2015l for a recent discussion), it seems 
to be due, mostly, to the effect of minor and major merg¬ 
ers teurke & i Gollina 20131), with minor mergers domin ating 
at later times ( Shankar et alJ l2013l : IBurke et alJ l2015l l . At 
most, BCGs may have grown by a factor ~ 1.8 in stellar 
mass since 2 ~ 1, although this factor could have been as 
small as ~ 1.2 if about half of the accreted stellar mass from 


the merging companion s became part of the intra-cluster 
light teurke et alJ 120151) . This mass growth seems to have 
been faster in the pas t, when both minor and major mergers 
were more common teurke fc Colling 120131 ). but these au¬ 
thors also found that BGGs in similar mass clu sters can have 
very d ifferent merging histories. Furthermore, lAscaso et al.l 
(I 2 OIIII reported that BGGs have grown in size by a factor 
of ~ 2 over a similar period. Interestingly, the difference 
in mass between cDs and elliptical BCGs in similar envi¬ 
ronments is of the order of 40% (i.e., comparable with the 
measured mass growth), and we find that the difference in 
size is a factor of ~ 2 (again, compatible with the measured 
size growth), but with a very large scatter in both cases. 
Additionally, IZhao et al.l (I 2 OI 5 II found that, when it could 
be reliably measured, the fraction of the light (stellar mass) 
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contained in the cD envelopes is of the order of ~ 40-60%, 
with significant galaxy-to-galaxy variations^ It is therefore 
plausible that most present-day BCGs started their life as 
ellipticals, and they subsequently grew, in stellar mass and 
size, due to mergers to become cDs. In this process, the 
characteristic cD envelope developed. The large scatter in 
the stellar masses and sizes of the cDs is explained by their 
different merger histories. Furthermore, the growth of the 
BCGs in mass and size seems to be linked to the hierar¬ 
chical growth of the structures they inhabit: as the groups 
and clusters become denser and more massive, the BCGs at 
their centres also grew. 

By the present tim e, most BCGs see m to be well ad¬ 
vanced in this process. IZhao et al.l ll2015ll found that the 
majority (~ 57%) of the BCGs are cDs, 21% have inter¬ 
mediate cD/E or E/cD morphologies, while ellipticals are 
a minority 13%). The presence of intermediate mor¬ 
phological classes suggests that this process is still ongo¬ 
ing. Present-day elliptical BCGs may (or may not) develop 
cD-type envelopes in the future, depending on whether the 
current merger rate is sufficient. With the limited statistical 
evidence that we have, we can only speculate about the ori¬ 
gin of the few 7%) BCGs with spiral and SO morpholo¬ 
gies, but perhaps these are the ones which avoided major 
mergers in their past history and retained their disks. 

If the evolutionary framework we propose is correct, 
one would expect the morphological mix of BCGs to change 
with redshift: at earlier times, the fraction of elliptical BCGs 
should be higher than today, with cDs showing the oppo¬ 
site trend. We have visually exami ned the images of t he 13 
BCGs in the ESO Distant Survey (I White et al.ll2005h clus- 
ters and groups f or which deep HST images are available 
dDesai et al.ll200'7tl . and morphologically classified them fol¬ 
lowing the same criteria used for the low-redshift sample. 
The average redshift of these galaxies is z 0.6. Although 
cosmological surface-brightness and resolution effects would 
have to be properly accounted for in a more systematic 
study, we feel that these HST images have enough resolution 
and depth (4 orbit exposure) for this purpose. They compare 
favourably with the SDSS images of the lower-redshift galax¬ 
ies. Notwithstanding these possible caveats, we find that 4 
of the BCGs are ellipticals, 3 cDs, 4 E/cD or cD/E, one is a 
spiral, and one is a merger. Although the sample is pitifully 
small, the trend seems to go in the right direction: the frac¬ 
tion of ellipticals more than doubles when compared with 
the local sample, while the fraction of cDs halves. There is 
also a significant fraction of galaxies with intermediate mor¬ 
phologies, suggesting that the transformation process is also 
happening at these redshifts. Of course, with such small sam¬ 
ple, no firm conclusions can be obtained, but at least these 
findings are compatible with our hypothesis. A systematic 
study of a large, well-defined sample of BCGs with deep 
HST images, reaching z 1, would be required to obtain a 
definitive answer. 

Numerical simulations and semi-analytic models (see, 

® Note that the galaxies for which this fraction could be reliable 
measured are the ones whose profiles are better modelled using 
two-component Sersic+exponential profiles. Since these tend to 
be the ones with more prominent envelopes, the average fraction 
of light in cD envelopes is probably closer to ~ 40%, the bottom 
end of the measured range. 


e.g., 1 ^ Lucia fc Blaizotl l2007l and references therein) pro¬ 
vide a plausible inside-out scenario for the growth of BCGs 
which is broadly compatible with our findings. At early times 
{z 1-3), dissipative processes similar to the ones proposed 
for the formation of normal giant elliptical galaxies were 
responsible for the building of the BCGs’ inner (elliptical¬ 
like) stellar component, whose light profile can be well rep¬ 
resented by a Sersic model. Subsequently, as the structures 
around BCGs grew hierarchically, the mass and size of these 
galaxies continued to increase, mainly due to dissipationless 
(dry) mergers, and the cD envelopes were formed as a result. 
This picture is also largely consistent with other observa¬ 
tions. For example, dry merg ers have been direct ly observed 
in cluster environments fe.g.. lvan Dokkumll2005ll . and it has 
been suggested that the accreted stars could built up the 
extended stellar halo s observed in BCGs llAbadi et al.ll20()6l : 
iMurante et al.l[2007h . 


5 CONCLUSIONS 

Using a large well-defined sample of 425 nearby Brightest 
Cluste r Galaxies from the catalogue of ivon der Linden et al.l 
ll2007lj . we have carried out a study of the relationships be¬ 
tween their internal properties (stellar masses, structural pa¬ 
rameters, sizes and morphologies) and their environment. 
The stellar masses M, are based on the MPA-JHU SDSS 
DR7 measurements. The structural parameters ( effective ra¬ 
dius R e and Sersic-index n) were derived by IZhao et al.l 
1 I 2 OI 5 II using single Sersic pr ofile fits. The visu al morpholo¬ 
gies were also published by IZhao et al.l (l2015l l , who found 
that the majority (~ 57%) of the BCGs are cDs, ~ 13% 
are ellipticals, 21% belong to intermediate cD/E or E/cD 
classes, and 7% have disk morphologies, with spirals and 
SOs in similar proportions. We use two separate measure- 
ments of the enviro nment, the local environmental density 
l|Tempel et al.l[2012l l. and the global dark-matter halo virial 
mass M 200 derived from t he cluster velocity dispersions 
llvon der Linden et al. lEo^. Our main conclusions are: 

• The Sersic-index n does not correlate with the stellar 
mass M, or the environment of the galaxies. 

• The effective radius Re of the BCGs correlates with 
their stellar mass M*, but the scatter is large (~ 0.3 dex in 
effective radius at a given mass). This correlation does not 
depend significantly on the environment. 

• Almost all BCGs have larger Re than non-BCG early- 
type galaxies of similar M*. The median radius of the BCGs 
is about twice as large as that of non-BCG early types of 
similar masses. This difference is largely due to the cD galax¬ 
ies, which dominate the sample. Moreover, the scatter in the 
Mt,-Re relation is significantly larger for the BCGs than for 
the other early-type galaxies, suggesting a more complex 
formation history. 

• More massive BCGs tend to inhabit denser regions and 
more massive clusters, but M* correlates significantly more 
strongly with environmental density than with the cluster 
dark-matter halo mass M 200 . Indeed, the apparent correla¬ 
tion between M* and M 200 can be explained by the correla¬ 
tions between M 200 and M* with environmental density. 

• The median stellar mass of cD BCGs is 2.1 x IO'^^Mq, 
~ 50% larger than that of ellipticals (1.4 x IO^^Mq). BCGs 
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with disk morphologies have even smaller stellar masses (me¬ 
dian 1.0 X 1O^^M0). 

• cDs seem to prefer marginally denser regions (by ~ 20% 
on average) than elliptical BCGs. Disk galaxies tend to live 
in the regions with the smallest densities. Similarly, cDs ap¬ 
pear to be hosted by more massive clusters/groups than el¬ 
lipticals (by factor of ~ 1.7 in median M 2 oo). However, these 
differences are only significant at the 2-2 .4(t level. 

• cD and elliptical BCGs show parallel correlations be¬ 
tween their stellar masses and environmental densities: 
larger and more massive galaxies tend to prefer denser en¬ 
vironments. However, at a fixed environmental density, cDs 
are, on average, ~ 40% more massive than elliptical BCGs. 
Due to the correlation between Re and M*, cDs and ellip¬ 
ticals also exhibit positive and parallel correlations between 
their effective radii and the environmental density. cDs are, 
statistically, twice as large as elliptical BCGs at a given den¬ 
sity. Disk BCGs tend to be smaller and less massive. 

Our results, together with the hndings of previous ob¬ 
servational and theoretical studies, suggest an evolutionary 
link between elliptical and cD BCGs. BCGs have experi¬ 
enced a significant growth in mass and size in the last ~ 6- 
8 Gyrs, largely due to the effect of minor and major mergers. 
The mass growth seems to have been faster in the past, when 
both minor and major mergers were more common, with 
minor mergers probably playing a dominant role in recent 
times. The amount of growth in mass and size experienced 
by BCGs since 2 : ~ 1 is comparable to the difference in mass 
and size between cDs and elliptical BCGs in similar environ¬ 
ments. Additionally, the fraction of the light (stellar mass) 
contained in the cD envelopes is also comparable with the 
average stellar mass difference between cDs and ellipticals. 
We therefore suggest that most present-day BCGs started 
their life as ellipticals, and they subsequently grew in stel¬ 
lar mass and size, due to mergers, to become cDs. In this 
process, the characteristic cD envelope developed. The large 
scatter in the stellar masses and sizes of the cDs is explained 
by their different merger histories occurring at 2 < 1. Fur¬ 
thermore, the growth of the BCGs in mass and size seems to 
be linked to the hierarchical growth of the structures they 
inhabit: as the groups and clusters became denser and more 
massive, the BCGs at their centres also grew. 

This process is nearing completion by the present time, 
since the majority of the BCGs in the local Universe have cD 
morphology. However, the presence of intermediate morpho¬ 
logical classes (cD/E and E/cD) suggests that the growth 
and morphological transformation of some BCGs is still on¬ 
going. It is also possible that today’s elliptical BCGs may 
develop cD-type envelopes in the future, depending on the 
merger activity they may experience. We also speculate that 
the BCGs with spiral and SO morphologies represent the mi¬ 
nority of BCGs which avoided major mergers in the past, 
thus retaining their disks. 

This scenario is broadly compatible with hierarchical 
inside-out models for the formation and growth of BCGs. 
Early dissipative processes were responsible for the build¬ 
ing of the BCGs’ inner elliptical-like stellar component. As 
the structures around BCGs grew hierarchically, the mass 
and size of these galaxies continued to increase, mainly due 
to dissipationless mergers, and the cD envelopes were thus 
formed. 


The evolutionary framework we propose seems to be 
able to explain the observed properties of BCGs, including 
the differences between the morphological classes. The obvi¬ 
ous next step to test this scenario is to carry out a study of 
the morphology, mass, structure and environment for a large 
and statistically robust sample of BCGs as a function of red- 
shift, reaching 2 ~ 1. A key piece of evidence would be the 
evolution of the fraction of cD BCGs with time, and its links 
with the growth of their masses, sizes and environments. 
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APPENDIX A: COMPARISON WITH GUO+OO 
A1 Structural parameters 

There are 104 galaxies in common between our sample and 
that of iGuo et alJ ll2009l l. A comparison between the mea¬ 
surements of the effective radius Re and the Sersic index 




log K,guo / kpc "Guo 

Figure Al. Comparison between the va lues of the e ffectiv e ra- 
dius Re and Ser sic index n obtained by IZhao et alJ ll2015^ and 
I Guo et al.l ll2009h for the 104 galaxies in common. The solid lines 
correspond to the 1-to-l relation. 

n for these galaxies is presented in Fig. lAll Although the 
measurements correlate very well, there are some relatively 
small systematic differences. T he median offset b etween our 
Re measurements and those of iGuo et al.l ll2009ll is 0.15dex. 
The median offset in n is 0.47. The larger values we obtain 
ar e due to imp r ovem ents in the sky subtraction implemented 
bv IZhao et akl (l2015lL In that paper we showed that the sky 
values provided by SDSS DR7 were overestimated due to the 
presence of extended objects. This is particularly important 
in crowded fields such as th e centres of groups and clus¬ 
ters. We used GALAPAGOS llBarden et al l2012l 'l to obtain 
a more reliable estimate of the sky after removing contam¬ 
ination from neighbouring objects. Although the reduction 
in the sky values is quite small (typically ~ 0.4 counts, or 
0.3%), the effect on Re and n can be significant for extended 
objec ts such as BGGs. M ore details are provided in section 
3.2 of IZhao et akl (l2015l L 

A2 Stellar masses 

In Section we found no correlation between n and M* 
fo r the BGGs i n our sample. This contrasts with the findings 
of IGuo et akl ll2009ll . who show a clear positive correlation 
in the sense that more massive BGGs seem to have higher 
values of n. In this Appe ndix we exp l ore th e possibility that 
the correlation found by I Guo et ^ (l2009l l may be due to 
the fact that these authors estimated stellar masses from the 
total luminosity derived from single Sersic model fits. These 
luminosities (and the derived stellar masses) are therefore 
model dependent, and, in particular, they will depend on 
the value of n. Since there is a direct relation between the 
best-fit t otal flux and n fo r a Sersic profile (see Equations 4 
and 6 in IPeng et al.l[201^ . this dependency could drive the 
observed correlation. 

In order to confirm this, we have derived stellar masses 
for the BGGs in our sample following the same method as 
IGuo et ^ (l2009l j using our own single Sersic fits. Since we 
have 104 BGGs in common with lGuo et akl ll2009lj . we can 
check that the values of M* derived in this way for the galax¬ 
ies in common agree well with theirs: the scatter in this com¬ 
parison is below 0.1 dex and there is no bias. In Fig. IA2I we 
show that, using these model-dependent M* values, a posi¬ 
tive correlation between n and M* is indeed found (Pearson 
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Figure A2. Sersic index n vs. stellar mass for the BCGs in our 
sample, similar to Fig.^ but with the stellar mass is derived 
following the method described in I Guo et al.l ll2009lb Symbols as 
in Fig. ^ See text for details. 


correlation coefficient 0.38). The correlation we fi nd is qual¬ 
itative ly similar to the one shown in Fig. 6 of I Guo et al.l 
(|2QQ9[) when considering the same mass range. 

This indica tes that the correlation claimed by 
IGuo et alJ (l2009l) may be the consequence of assuming that 
a Sersic model fit provides an accurate representation of 
the total light distribution of BCGs. This assumption is 
clearly not correct, particularly for cD galaxies, as demon¬ 
strated by previous studies (see lZhao et al.l[2015l and refer¬ 
ences therein). Measuring the total luminosity of a galaxy 
is far from trivial and, of course, the Petrosian magnitudes 
used to derive MPA-JHU masses are not without their prob¬ 
lems (see, e.g.. iGraham et 'aDbOOsI l. We do not claim that 
the stellar masse s we use are better than the ones used by 
I Guo et ^ ll2009ll . but they are, at least, model independent 
and not directly linked to the models used to derive the 
structural parameters that we study. For these reasons we 
prefer to use the MPA-JHU masses in this paper. Neverthe¬ 
less, bearing in mind this uncertainty, we have checked and 
confirmed that all our conclusions (with the exception of the 
lack of correlation between M* and n) remain the same if we 
use Sersic-model based luminosities/stellar masses instead of 
the MPA-JHU ones. 





























